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Hybrid lateral superlattices composed of a square array of antidots and a periodic one-dimensional 
magnetic modulation are prepared in Ga[Al]As heterostructures. The two-dimensional electron 
gases exposed to these superlattices are characterized by magnetotransport experiments in van¬ 
ishing average perpendicular magnetic fields. Despite the absence of closed orbits, the diagonal 
magnetoresistivity in the direction perpendicular to the magnetic modulation shows pronounced 
classical resonances. They are located at magnetic fields where snake trajectories exist which are 
quasi-commensurate with the antidot lattice. The diagonal magnetoresistivity in the direction of 
the magnetic modulation increases sharply above a threshold magnetic field and shows no fine struc¬ 
ture. The experimental results are interpreted with the help of numerical simulations based on the 
semiclassical Kubo model. 

PACS numbers: 73.23.-b, 73.63.-b 


I. INTRODUCTION 

Artificial lateral superlattices (LSLs) in two- 
dimensional electron gases (2DEGs)i - — are of great 
interest for fundamental studies of the electron dynamics 
in periodic potentials. Since it is very common that 
the artificial lattice constants place the systems in the 
transition region between the quantum and the classical 
regime, classical, semiclassical as well as quantum 
descriptions are all justifiable and enable studies of 
the validity of these approaches including their limits. 
Besides the Fermi wavelength \p and the electronic co¬ 
herence length, the elastic mean free path is an important 
parameter as well, since it defines the length scale below 
which interaction with the LSL potential dominates 
over random scattering. Many different variants of LSLs 
have been investigated in great depth. One-dimensional 
(Id) electrostatic)^ and magnet io^ - — lattices, where 
the modulation extends along one spatial coordinate 
and the structure is homogeneous along the second 
coordinate, show magnetoresistivity resonances that can 
be explained in terms of guiding center drift resonances 
of the cyclotron motion within a classical picture^ or 
by miniband formation in a quantum picture ^ 10 ! 11 
One-dimensional magneto-electric hybrid LSLs have 
been studied in some experiments as well, where the 
strain imposed by the ferromagnetic or superconductive 
electrodes used to define the magnetic LSL also gen¬ 
erates an electrostatic super latticed Two-dimensional 
LSLs, both magnetioiii^ and electrostatic^ - — , have 
been studied in thoroughly as well. Their classical 
dynamics corresponds to a mixed phase space where 
chaotic and regular dynamics coexis t 14 ! 15 and causes 


commensurability resonances that are characteristic for 
the type of Bravais super lattice employed, like square^ - —, 
rectangular— or hexagona l 17 ! 18 . Within a quantum 
picture, on the other hand, a fractal energy spectrum, 
also known as Hofstadter butterfly ^ 19 ! 20 is seen for weak 
electrostatic modulation amplitudes^ B-periodic os¬ 
cillations on top of commensurability resonance s 16 ! 18 ! 22 
can be explained within a semiclassical approach by the 
Aharonov-Bohm22 - or Altshuler-Aronov-Spivab^i- effect 
in terms of quantized motion along closed trajectories 
defined by the LSL potential and the magnetic fields. 

2DEGs with very high electron mobilitie s 26 ! 27 have re¬ 
cently been developed into mature systems. They enable 
the preparation of LSLs with large lattice constants in 
the classical ballistic regime and facilitate the definition 
of novel types of LSLs with more complex unit cells. 
Here, the study of such a hybrid LSL, composed of a 
a two-dimensional, square antidot lattice and a one¬ 
dimensional magnetic array is reported. The magnetic 
LSL consists of approximately Lorentzian shaped peaks 
of alternating sign and thus has a vanishing average 
magnetic field. Snake trajectories, i.e., trajectories 
formed by the superposition of an oscillatory motion 
along the first direction and a motion with nonvan¬ 
ishing average velocity along the second direction^ 
can become commensurate with the antidot lattice, 
and magnetoresistivity resonances are to be expected. 
Furthermore, for the magnetic modulation amplitudes 
applied here, closed electronic orbits are absent. 

After the sample preparation and the experimental 
setup are introduced in Section mi the measurements 
are presented in Section [TTT] and interpreted with the 
help of numerical simulations in Section [iVl The paper 
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FIG. 1. (color online), (a) Scheme of the hybrid lateral su¬ 
perlattice geometry. Circular holes (red circles) that form the 
antidots are etched into the sample. Every second line of an- 
tidots is covered by a Dy stripe of 1 pm width. The top view 
of the pattern is shown in the uppermost part. In the mid¬ 
dle, a cross section in the yz - plane at x=0 is shown. The 
2DEG is depleted underneath the etched regions. The Dy 
stripes are magnetized in the y - direction, as indicated by 
the arrows. The corresponding electrostatic potential (possi¬ 
ble strain effects are neglected) and the perpendicular mag¬ 
netic field B z (y ) are shown in the lowermost part, (b) Sketch 
of the sample layout. Two identical hybrid LSLs, (i) and (ii), 
are defined in an L-shaped Hall bar. The two components of 
the hybrid LSL, namely the antidot lattice (iii) and ferromag¬ 
netic stripes (iv) are defined separately, and the edge of a Dy 
pad centered in a Hall cross (v) enables Hall magnetometry. 
The coordinate system shows the crystallographic orientation 
of the Hall bar. (c) p xx of array (iii) as a function of a homo¬ 
geneous perpendicular magnetic field, (d) p xx (B y ) of array 
(iv). 


concludes with a summary and an outlook in Section [VI 


II. SAMPLE PREPARATION AND 
EXPERIMENTAL SETUP 


A GaAs/Alo. 3 Gao. 7 As heterostructure with a 2DEG 
90 nm below the surface is used. After a brief illumi- 
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where zq is the distance between the 2DEG and 
the bottom of the Dy film, j is an integer and N 


nation with infrared light, the unpatterned 2DEG has a 
density of 3.6 x 10 15 m -2 and a mean free path of 88 p m 
at liquid helium temperatures. The sample geometry 
is sketched in Figs, [lja) and (b). An L-shaped Hall 
bar, oriented parallel to the natural GaAs cleavage direc¬ 
tions, was prepared by optical lithography. Three identi¬ 
cal, square antidot lattices (lattice constant a = 1.0 pm) 
were patterned on one Hall bar by electron beam lithog¬ 
raphy and subsequent reactive ion etching. Lithographic 
antidot diameters of duty = 200 nm (sample A) as well 
as duth = 300 nm (sample B) were prepared on sepa¬ 
rate Hall bars. As a consequence of a lateral depletion 
length of 45 nm around the antidots, this corresponds 
to electronic diameters of d ~ 290 nm and d ~ 390 nm, 
respectively, as measured by the Aharonov-Bohm oscil¬ 
lation period observed in large magnetic fields j 29 i 30 Since 
(a — d) / Xf ~ 17 fro sample A and ~ 14 for sample B, 
respectively, these LSLs reside well inside the classical 
regime. After the definition of the antidots, Dy stripes 
of width a and a period of 2 a were prepared on top of 
two antidot lattices by electron beam lithography, en¬ 
abling measurements of all resistivity components in one 
cooldown, see Fig. Hlb). The Dy stripes have a thickness 
of h = 250 nm to ensure a strong fringe field when mag¬ 
netized. In sample A, they were deposited directly on the 
GaAs, while in sample B, a homogeneous film of 5 nm Cr 
plus 5 nm Au thickness was evaporated on top of the an¬ 
tidot lattice prior to the Dy deposition. This allows us 
to estimate the role of strain effects^^i possibly induced 
by the Dy stripes, which are centered at the columns of 
antidots and aligned parallel to the x-direction. The lat¬ 
eral size of the super lattices is 100 pm in longitudinal and 
50 pm in transverse direction (100 x 25 unit cells). For 
control measurements, the Hall bar furthermore contains 
a nominally identical magnetic stripe array without the 
antidots underneath, and the edge of a Dy film in a Hall 
cross for Hall magnetometry j 32 1 33 

The samples were inserted in a 4 He gas flow cryostat with 
a variable temperature insert and a base temperature of 
1.4 K. The system is equipped with a magnet of 8 T max¬ 
imum field strength. The external magnetic field B y was 
applied in the y - direction. It magnetizes the Dy stripes 
to a magnetization of poM(B y ). The 2DEG responds 
predominantly to the z-component of the fringe field of 
the Dy stripes, and we therefore neglect the influence of 
in-plane magnetic fields on the 2DEG throughout this 
paper. The magnetic field profile B z (y) is indicated in 
the lowermost section of Fig. HJa). From the fringe field 
of a perfectly magnetized stripe, one expects^ 


[y ~ a (2j T ^)] 2 + 

[y - a(2j T \)] 2 + (z 0 + h) 2 


( 1 ) 


denotes the total number of Dy stripes. This magnetic 












































3 


profile has peaks of alternating sign with amplitude 
B™ ax {B y ) = \B z {y=[2j-\]a,B y )\, The maxi- 

mum magnetization of our Dy films is /io M ~ 2.7 T 
for By > 5 T, corresponding to an upper limit of 
-grriax ^ 480 mT. The coercive magnetic field is 
B c = 670 mT. The resistivity components pij(B) with 
i,je{x,y} were determined by applying an AC current 
of 100 nA with a frequency of 17.7 Hz from source 
S to drain D, see Fig. mb), and by measuring the 
electrostatic potentials in x- and y-direction at voltage 
probes with a lock-in amplifier. 


III. EXPERIMENTAL RESULTS 




FIG. 2. (color online), (a) Magnetoresistivities p xx (B y — B c ) 
of the hybrid superlattices of samples A and B, shown for the 
up-sweeps (increasing B y ). The inset shows the hysteretic 
behavior, exemplified for sample A. The dashed lines indicate 
B c . (b) Temperature dependence of p xx (B y ) for sample A, 
as observed in a different cooldown. Only the down-sweep 
direction is shown for clarity. The measurement of p yy (B y ) 
for sample A is reproduced in (c). 

In Figure [l](c) and (d), the magnetoresistivities of the 
LSL components of sample A for the two hybrid LSL 
components, namely the square antidot lattice (iii) and 
the array of magnetic stripes (iv), respectively, are re¬ 
produced. The antidot lattice reveals the well-known 


commensurability resonances with resistivity maxima at 
perpendicular magnetic fields where the cyclotron or¬ 
bit is commensurate with one, two, four or nine en¬ 
closed antidots^ - — For B z > 250 mT, Shubnikov - de 
Haas oscillations set in. p X x{B y ) of the Dy stripes 
shows a peak centered at B c and some weak features 
at larger magnetizations. This type of magnetoresis¬ 
tivity of magnetic stripe arrays in in-plane magnetic 
fields has been studied theoretically^ while to the best 
of our knowledge, experiments have been reported only 
in related configurations^ Numerical simulations based 
on the classical Kubo formalism (see below for details) 
give a weak, positive magnetoresistivity without fine 
structure^ as measured for \B y \ > 4T. This indicates 
that the peak at B c is not an intrinsic classical prop¬ 
erty of the magnetic profile itself, and we tentatively at¬ 
tribute it to the frequently observed and still not fully 
understood negative colossal magnetoresistance in high- 
mobility 2 DEGsj 2I— 9 which is beyond our focus here, 
possibly in combination with other effects like weak lo¬ 
calization. The strength of this feature depends on the 
cooldown cycle. It should be emphasized that p xx of the 
2DEG underneath the Dy array is constant over the full 
scan range within ±0.8 D. For the following, this contri¬ 
bution can therefore be neglected. 

The diagonal magnetoresistivities p xx {B y ) and 
p yy (B y ) of the hybrid LSLs are reproduced in Fig. [2j 
and we first focus on p xx {B y ) as observed on samples A 
and B (a). As B y is detuned away from F> c , a positive 
magnetoresistivity is observed. As B y is further in¬ 
creased, two peaks are seen, separated by a pronounced 
minimum. Around B y — B c ^ 600 mT, a decrease 
of p xx by roughly a factor of 2 is seen, followed by a 
broad maximum that extends up to B y — B c « 2T. 
These most prominent features appear in sample B 
at somewhat larger magnetic fields than in sample A. 
Also, even though the positive magnetoresistance is less 
pronounced in sample A than in sample B, sample A 
shows clear additional finer structures, some of which 
are also adumbrated in p xx (B y ) of sample B. These 
differences can be traced back to the Cr/Au electrode 
present in sample B, as will be discussed below in 
more detail. In the following, we focus on sample A. 
In the inset of Fig. [2] (a), the hysteretic behavior of 
Pxx{B y ) is reproduced. The features are fairly symmetric 
about H c , while the symmetry of the up-sweep to the 
down-sweep about B y = 0 is close to perfect. This 
behavior indicates that the magnetization of the Dy 
stripes is not perfectly antisymmetric about B y = B c 
(see below). These magnetoresistivity features show a 
weak temperature dependence, see Fig. [2] (b), and the 
most pronounced ones remain visible up to T « 16 K. 
This suggests that they should be interpretable within 
a classical picture. They are furthermore superimposed 
to a slowly varying negative magnetoresistivity that 
extends to \B y — B c \ « 1.6T, after which it increases 
slightly. This background depends somewhat on the 
cooldown cycle. The strong positive magnetoresistivity 
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in a narrow interval around B c is still clearly visible 
at 32 K, and behaves similarly to that one observed in 
two-dimensional antidot lattices, see also Fig. Q](c). It 
is due to a B z - induced increase in scattering at the 
antidots and is of no further interest here. 

A smooth increase of p yy (B y ) is observed as B y is 
driven away from B c . A sharp increase sets in for 
|By — B c | « 1.6T and stops for | B y — B c \ « 3.5T, see 
Fig. [2] (c). The shape of p yy (B y ) strongly resembles 
that one observed for single magnetic barriers^ as well 
as magnetic barriers in series of alternating polarity4£ 
Within a classical picture, the increasing amplitude of 
B z (y ) reflects an increasing fraction of the incident elec¬ 
trons that gets reflected at the magnetic barrier. Above 
a threshold amplitude of B z (y ), electrons can only pass 
the barrier via Ex B drift at the edges of the Hall bar, or 
by scattering events inside the magnetic barrierThese 
effects cause a saturation of p yy at large B™ ax . Since 
our ferromagnetic array represents an array of magnetic 
double barriers in series^, p yy (B y ) can thus be qualita¬ 
tively understood in terms of the properties of magnetic 
double barriers with the antidots acting as scatterers^S 
and is not a unique signature of the the hybrid lattice. 
The onset of the sharp increase of p yy furthermore cor¬ 
relates with the end of the negative magnetoresistivity 
in x-direction. Comparison of p xx {B y ) to Pyy{B y ) re¬ 
veals that the transport at large magnetic fields is highly 
anisotropic. For example, for | B y — B c \ = 2.5 T, the ra¬ 
tio p yy /pxx reaches a value of « 230. This suggests that 
for sufficiently large magnetization of the Dy stripes, the 
electrons are guided along the x-direction by the mag¬ 
netic modulation, while crossing the magnetic walls is 
highly unlikely. 

The off-diagonal components of the magnetoresistivity 
tensor were measured as well4£ Since the average per¬ 
pendicular magnetic field is zero, they vanish to a good 
approximation in the magnetic field range where the res¬ 
onances in p xx (B y ) appear and are thus not very helpful 
for their interpretation. 


IV. MODEL CALCULATION AND DISCUSSION 

A coarse estimation, to be substantiated below, reveals 
that for B™ ax < 500 mT, B z (y) is too weak to generate 
closed cyclotron-type orbits. Therefore, the magnetore¬ 
sistivity resonances must originate from open trajecto¬ 
ries. This situation is quite different in comparison to 
antidot lattices in homogeneous magnetic fields where 
closed orbits, runaway trajectories and chaotic orbits co¬ 
exist and all contribute to the magnetoresistivity with a 
magnetic field-dependent weight^ Open cycloid orbits 
are absent as well in the interval where the resonances ap¬ 
pear, and it is therefore expected that snake trajectories 
play an important role, the most obvious type of which 
is centered at the roots of B z (y) and runs along columns 
of antidots in x-direction. Since this is a classical picture 


and moreover the most pronounced features of the mag- 
netoresistivity show a weak temperature dependence, it 
appears plausible to model them using the classical Kubo 
formalism. The code we use has been presented in detail 
elsewhere^ and is therefore only briefly sketched here. 
We show the simulations for the parameters of sample 
A. Electrons are injected at random locations inside a 
unit cell of the LSL. They initially move in random di¬ 
rections with their Fermi velocity of vp = 2.6 x 10 5 m/s. 
The incremental change of the direction of motion by the 
inhomogeneous magnetic field given by Eq. 0 is calcu¬ 
lated with a step width of 2 nm, and specular reflection 
at the antidots with d = 290 nm is assumed. Further¬ 
more, zo = 90 nm is used, and we assume that the an¬ 
tidot potential is hard-wall, as justified by the large a/d 
ratio. The simulations are carried out for zero tempera¬ 
ture. From the simulated diffusion tensor obtained via 
the Kubo formula, the magnetoresistivity components 
are obtained via the Einstein relation for a degenerate 
2DEG. 

Figure [3] (a) shows the simulated magnetoresistivity 
p xx as a function of the maximum of the perpendicular 
magnetic field B ™ ax , see also Figs. [I] (a) and (b). 



FIG. 3. (color online), (a) p xx (B™ ax ) as simulated within 
the Kubo model, plotted as a function of the maximum of 
B z . The period of the commensurate snake trajectories is 
indicated at the most prominent resistivity minima. The ex¬ 
perimental trace of sample A has been scaled to B z with the 
help of the magnetization trace, as obtained from the Hall 
resistance of one Dy edge centered inside a Hall cross (b). (c) 
Simulation results for p yy (B™ ax ) in comparison to the scaled 
experimental data of sample A. 
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As in the experiment, several features in p xx are ob¬ 
served. Close to B™ ax < 250 mT, a positive magnetore¬ 
sistivity is present. For B™ ax < 250 mT, a series of re¬ 
sistivity minima at B™ ax = 32 mT, 53 mT and ~ 110 mT 
is visible. A clear but weaker additional minimum is vis¬ 
ible at B™ ax « 170 mT. Above a sharp decrease of p xx 
at B™ ax « 260 mT, a broad minimum around 280 mT 
is present, followed by some weakly pronounced max¬ 
ima and minima. Finally, another sharp decrease of p xx 
around B™ ax = 500 mT is observed. A direct compar¬ 
ison with the measurements requires knowledge of the 
transformation function poM(B y ). Conceptually, it can 
be determined by Hall magnetometry of the stripe array 
on top of Hall crosses well inside the diffusive regime. 
In the ballistic or quasi-ballistic regime, the Hall volt¬ 
age translates into the magnetization by nontrivial cor¬ 
rection factors ^ 33 ! 44 the detailed discussion of which is 
beyond our scope here. Since such an estimation would 
still assume perfect, mono-domain magnetization of the 
Dy stripes as well as a certain shape of the fringe field, 
some uncertainty would remain. Therefore, in order to 
estimate poM(B y ), we restricted ourselves to Hall mag¬ 
netometry of the edge of a Dy film, prepared in the same 
process step as the magnetic lattice. The measured Hall 
voltage as a function of B y , reproduced in Fig. [3] (b), 
shows a marked peak where the average cyclotron diame¬ 
ter equals the width of the voltage probe. The decrease of 
the Hall voltage at larger magnetic fields originates from 
ballistic effects^ The asymmetry of the Hall voltage fur¬ 
thermore indicates that the magnetization of the film is 
not perfect. Therefore, we compare the measured data 
to the simulations by scaling it with an approximated 
function poM(B y ), obtained numerically along the lines 
of Ref. |33, where p^M is roughly proportional to B y for 
| B y — B c | < 600 mT and depends only weakly on B y for 
larger applied magnetic fields. This analysis of the Hall 
magnetometry indicates a saturation magnetization for 
the Dy stripes of « 2.7 T, and B c = 670 mT can be read 
out directly. The data measured at sample A in the up¬ 
sweep for B y > B c in Fig. [2] (a) are scaled accordingly 
and replotted in Fig. [3] (a) as a function of B ™ ax , which 
allows a more direct comparison to the simulations. 
Even though the simulated function p xx (B™ ax ) devi¬ 
ates from the experimental trace in several aspects, the 
most prominent features are reproduced qualitatively, 
namely the positive magnetoresistivity around B™ ax = 
0, minima close to B™ ax = 53 mT, 110 mT, 170 mT and 
280 mT, the decrease of p xx at B™ ax « 260 mT, and 
some weakly pronounced maxima and minima at larger 
magnetic fields. The sharp decrease of of p xx around 
j^rnax _ 500 no t observed experimentally, most 

likely because our fringe fields are too weak. 

The simulation of p yy (B ™ ax ) is compared to the scaled 
experimental data in Fig. [3](c). Very good agreement is 
found for B™ ax < 0.3 T, while the strong increase of the 
resistivity around B™ ax « 0.4 T is reproduced as well, 
though shifted to slightly higher magnetic fields. Further 
simulations^ show that the presence of the antidots 


does influence p yy (B ™ ax ) to some extent, but the overall 
behavior is dominated by the magnetic barriers and is 
not an effect of the hybrid superlattice. 

We proceed by interpreting the magnetoresistivity fea¬ 
tures in terms of the electron dynamics which determines 
the components of the magnetoconductivity tensor 
The off-diagonal elements a xy and a yx are approximately 
independent of B™ ax and of the order of 0.1 mS. a yy 
decreases from 18 mS at B™ ax = 0 to almost zero at 
B™ax ~ 0.5 T. Only a xx shows resonances as B™ ax 
is changed. This implies that p xx ~ 1 /a xx and p yy ~ 

l/CTyy, While p X y(By) & (T~y / {(T XX (Tyy) . TllC Sl^Up Ul- 

crease of p xy (see Ref. 0) and p yy at B™ ax = 0.5 T 
has thus its origin in the strongly suppressed diagonal 
conductivity in y-direction. 



FIG. 4. (color online). Poincare sections for various values of 
B™ ax (a-e). Some characteristic trajectories are shown in (f), 
the initial conditions of which are indicated by full circles in 
the corresponding Poincare sections. 

A deeper insight into the underlying electron dynamics 
can be gained by looking at characteristic electron trajec¬ 
tories. They can be identified with the help of Poincare 
sections, which illustrate the dynamics of the electrons 
by their coordinates in a (i, Vj) cross section of the phase 
space (i = x,y and Vj , j = x,y denote the position and 
velocity coordinates, respectively). We start with a dis¬ 
cussion of the minima of p xx at smaller magnetic fields, 
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B™ax < 280 mT. Each dot in Figs. 0](a-e) represents the 
coordinates of an electron passing with v x > 0 through 
one of the (y,v y ) - planes located at at x = raa, where m 
is an integer. The Poincare section for B™ ax = 265 mT 
(a) shows a pronounced accumulation of the electrons in a 
semicircle-like structure that extends over 85% of possible 
v y components. This region hosts quasi-commensurate 
snake trajectories with a wavelength very close to 2a. 
They run parallel to the magnetic stripes, as illustrated 
by the sample trajectories shown Fig. IU(f), and typically 
get scattered at the antidots after less than 30 snake pe¬ 
riods. Likewise, the Poincare sections for the minima of 
p xx at B™ ax = 170mT (b), 108mT (c), and at 53mT 
(d) reveal that here, quasi-commensurate snake trajecto¬ 
ries of various periodicity exist. They extend along the 
x - direction, and their weight decreases as the magnetic 
field is decreased, which correlates with the magnitude of 
the corresponding resistivity dips. Outside the resistiv¬ 
ity minima, such an accumulation of electrons in snake 
trajectories is not seen in the Poincare sections, as illus¬ 
trated for B™ ax = 240 mT in Fig. 0] (e). 

In addition, snake orbits exist which run at an angle % 0 
to the x-direction, as exemplified in Fig. 0] (f). In the 
Poincare sections, such trajectories form white regions, 
since the electrons do not return to the column in which 
the electrons start. They can be found over the whole 
interval where resonances are observed, and we do not 
find a correlation between their weight in the Poincare 
section and the magnetoresistivity. We furthermore em¬ 
phasize that, as anticipated above, closed orbits are ab¬ 
sent. It thus emerges that the minima of p xx {B y ) for 
B™ ax < 280 mT correlate with the presence of quasi- 
commensurate snake trajectories that run parallel to the 
magnetic stripes for many antidot periods, while snake 
trajectories running in other directions do not show such 
a correlation. Both the depletion and accumulation re¬ 
gions of the Poincare sections are embedded in an ap¬ 
proximately homogeneously filled background, which is 
due to electrons that move in snake orbits as well, but 
experience frequent scattering at the antidots. Typically, 
such trajectories complete no more that two snake peri¬ 
ods before they get scattered* 30 We note that both the 
accumulation and the depletion regions contain mostly 
not perfectly periodic trajectories and are thus chaotic 
as well. Regular orbits should exist inside the accumu¬ 
lation regions, but we have been unable to identify such 
points in the Poincare sections, which indicates that the 
regular regions have a very small volume. The composi¬ 
tion of the phase space of this hybrid LSL is thus different 
to that one of antidot lattices where disjunct, extended 
regular and chaotic regions coexist. 

It is remarkable that adjacent resistivity minima some¬ 
times correlate with accumulations of snake orbits of the 
same periodicity. For example, the minima at B™ ax = 
170 mT and at 108 mT both correlate with the accumula¬ 
tion of snake trajectories with a period close to 4a. While 
the snake trajectories that belong to the pronounced min¬ 
imum of p xx at 108 mT remain commensurate over a rel¬ 


atively large interval of magnetic fields and initial condi¬ 
tions, those found at the weak minimum at 170 mT, like 
the two shown in (f) with their location indicated in the 
Poincare section in (b), are more fragile. 

For B™ ax > 280 mT, the simulation shows a series of 
weakly pronounced features that end with a strong de¬ 
crease of p xx at B™ ax « 500 mT. Qualitatively similar 
features are observed experimentally for sample A and 
can be only guessed for sample B. In this interval, the 
Poincare sections show a rich pattern of accumulation 
regions, together with a few depleted areas, see Fig. [5] 
(a). This pattern evolves smoothly as a function of B™ ax 
without changing its qualitative appearance. The snake 
trajectories in this interval have a period close to 2a and 
show only a few oscillations before they get scattered. 
We also find occasional trajectory segments of skipping 
orbits, see Fig. [5] (b). Thus, the magnetoresistivity fea¬ 
tures in the interval 280 mT < B™ ax < 500 mT do not 
correlate in a straightforward way with characteristic tra¬ 
jectories. 



FIG. 5. (color online), (a) Poincare section for B™ ax = 
375 mT. (b) Some typical trajectories at B™ ax — 375 mT 
with initial conditions as indicated by full circles in (a), and 
a cycloid trajectory for B™ ax = 650 mT. 

The limit of large magnetic fields is characterized by 
B™ ax ~ 500 mT. It is not experimentally accessible in 
our samples. The simulations suggest that the decrease of 
p xx originates from the formation of cycloid trajectories 
which drift along the magnetic field peaks. An example 
of such a trajectory is shown in Fig. [5] (b). For B™ ax > 
520 mT, cycloid orbits exist that never hit an antidot. 
Therefore, highly conductive channels in x-direction are 
formed. 

We conclude this section by discussing possible reasons 
for the differences observed between sample A and sam¬ 
ple B, as well as for the deviations between the simula¬ 
tions and the experiments. The most prominent features 
in sample B appear at higher magnetic fields than their 
respective counterparts in sample A. With the support 
of the corresponding numerical simulations, this can be 
traced back to the larger distance of the Dy stripes to the 
2DEG due to the Cr/Au layer in between, which makes 
a higher magnetization necessary to achieve a fringe field 
of the same magnitude at the depth of the 2DEG. Also, 
the less prominent features observed at sample A are sup¬ 
pressed in sample B. This may be due to the larger an- 
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tidot diameters in sample B which is known to smear 
out commensurability resonances 4. Also, gating of high 
mobility heterostructures can decrease the mobility^ 
The simulated amplitudes of the commensurability 
oscillations are furthermore stronger than the measured 
ones. We attribute this partly to the deviations of the 
real magnetic field profile B z (y ) from the simulated 
one, which is to be expected from the asymmetric 
magnetization characteristics of the Dy film. Deviations 
from the assumed hard-wall potential may deform the 
trajectories, thereby weakening the resonances. Another 
possible reason are piezoelectric effects due to strain 
imposed by the Dy stripes, which could modulate the 
electron density and the mobility for our crystalli- 
graphic orientation of the Hall bars. This effect has 
been reported in the literature to get attenuated by 
depositing the stripes on top of a homogeneous metallic 
layer. 7 Therefore, by comparing the measurements of 
sample A with those of sample B, we conclude that 
if strain effects were relevant, they would generate 
additional fine structure rather than smearing them 
out. To further elucidate this issue, we have performed 
numerical simulations as described above, with an 
additional electrostatic potential of a cosine shape in y 
direction with the period of the magnetic stripes and 
a rather strong amplitude of 1 meV, in accordance to 
typical values found in the literature] 31 ! 45 Somewhat 
surprisingly, we do not find significant deviations of 
the resistivity from the unmodulated case (not shown) 
and therefore conclude that the magnetic field gradient 
dominates over electrostatic effects in the regime where 
the resonances are observed. Strain effects thus do not 
play a prominent role. Also, the simulation neglects 
finite size effects. For example, a magnetic barrier close 
to the Hall bar edges induces E x B drift, and electron 
scattering at the Hall bar edges may provide additional 
conductance channels. Finally, at the large in-plane 
magnetic fields present in our implementation, magnetic 
mass effects can deform the snake trajectories to a small 


extent^ 


V. SUMMARY AND CONCLUSIONS 

Hybrid magneto-electric lateral superlattices com¬ 
posed of a two-dimensional antidot array and a one¬ 
dimensional magnetic modulation have been defined in 
high-mobility two-dimensional electron gases and studied 
by transport experiments in a configuration with vanish¬ 
ing average perpendicular magnetic field. Despite the 
absence of closed trajectories, pronounced classical mag- 
netoresistivity resonances have been observed. The mag¬ 
netoresistivity minima correlate with the accumulation 
of electrons in snake trajectories, as observed in Poincare 
sections, that are quasi-commensurate with the antidot 
lattice and oriented along the direction in which the mag¬ 
netic field is homogeneous. Snake trajectories running in 
other directions are present as well, but their appear¬ 
ance does not correlate with the resistivity minima. The 
Poincare sections do not show extended regular islands. 
We hope that these rsults will trigger quantum simula¬ 
tions of this system which should be able to interpret 
the magnetoresistivity resonances on a more fundamen¬ 
tal level. The longitudinal magnetoresistivity is further¬ 
more strongly anisotropic, with resistivity ratios above 
200 for large magnetic fields. To a good approximation, 
however, the magnetoresistance in the direction perpen¬ 
dicular to the magnetic stripes can be understood as a 
resistance of magnetic barriers in series and does not re¬ 
veal super lattice-specific properties. Further experiments 
may comprise the application of additional homogeneous 
perpendicular magnetic fields, a more detailed study of 
p yy , the interaction of the electrons in snake trajecto¬ 
ries with resonant electromagnetic radiation, or magnetic 
mass effects. 
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